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TECHNICAL NOTE NO. 1026 

CHARTS OF THERMODYNAMIC PROPERTIES OF FLUIDS ENCOUNTERED IN 
CALCULATIONS OF INTERNAL COMBUSTION ENGINE CYCLES 
By* H. C. Hottel and G. C. Williams 

SUMMARY 


A single chart has teen prepared that simplifies the 
calculation of the thermodynamic properties of air, various 
octene-air mixtures, and mixtures of those with their prod- 
ucts of combustion at all temperatures "below which chemical 
dissociation becomes unimportant. The chart is based on the 
use of 1 pound mol of mixture, and examples of its use are 
given in the form of the calculation of a turbocompressor 
power plant and a supercharged Otto engine cycle. 


INTRODUCTION 


Reference 1 presented thermodynamic property charts for 
the products of combustion of octene and air in different 
ratios. For the calculation of various engine cycles, there 
is necessary in addition a chart or set qf charts giving the 
thermodynamic properties of unburned air- fuel mixtures. Such 
an unbur ned-fuel-air mixture chart has been calculated, and 
so plotted that it applies to various air-fuel ratios and 
various ratios of burned to unburned gases, and includes pure 
air as the basic chart , 

This investigation, conducted at the Massachusetts 
Institute of Technology, was sponsored by and conducted with 
the financial assistance of the National Advisory Committee 
for Aeronautics. 

The completion of this project would not have been pos- 
sible without the conscientious and efficient assistance of 
Mrs. Edward Addelson, who made all the calculations and did 
all the drafting and lettering on the final charts. 



FUNDAMENTAL DATA 


The charts in reference 1 were each on a basis of 1 
pound of air plus the appropriate fuel, and the coordinates 
were entropy S and internal energy E, A eet of correspond- 
ing "unburued" chart b on the same has! a could be prepared, 
one corresponding to each ‘’burned” chart. If, however, the 
basis for the properties of unburned mixtures is not 1 pound 
of air but 1 pound mol of mixture, the p-v-T relations for 
various mixture ratios are substantially the same - identical, 
of course, if the various air-fuel mixtures were perfect 
gases . 


Furthermore, if the coordinates of the chart are entropy 
S and temperature T rather than S and 3 as previously, 
the representation of the properties of different mixtures on 
a single chart is facilitated. It will be shown that the 
small variations, with air-fuel mixture ratios, from a common 
p-v-T relation due to gas imperfections and from a common H-T 
or E-T relation due to molal- speci f i c-heat variations from 
fuel to air, can be allowed for on a single chart. Such a 
chart can then be used for all mixtures, burned or unburned 
and including pure air, at the temperatures of interest here 
(up to 1800° R (1340° F) for unburned fuel-air mixtures, 

2520° R (2060° F) for burned mixtures). 

Although such a single chart has the disadvantage over 
a series of charts that a change in basis is involved during 
the calculation of some thermodynamic cycles, that change is 
not difficult (an example will be given;; furthermore, the 
advantages of a single chart are many. Besides the physical 
simplification of handling a single chart there is the possi- 
bility of calculation of complete gas turbine cycles of any 
mixture ratio of present interest instead of only those ratios 
for which burned charts are available; a comparison of the 
air-standard cycle with- the actual fuel cycle can be made di- 
rectly on the single chart; the chart gives directly the prop- 
erties of nure air for calculation of various pumping or 
power cycles. 

In the previous charts zero values were assigned to va- 
rious thermodynamic properties at a base of 100° F. Since 
calculations of high-altitude power-plant performance will 
involve the present chart at temperatures to -60° F or lower, 
it is desirable to drop the temperature base to a value 
below any anticipated in order to involve differences of pos- 
itive numbers in engineering calculations. The base chosen 
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was 200° K (360° R or -100° P) because basic thermodynamic 
data are tabulated in the literature at that temperature. 


DESCRIPTION OP CHART 
Used as an Air Chart 


The chart is best viewed first without consideration of 
the various peripheral plots. It is then a representation 
of the p-v-T-S relations for pure air, -due allowance having 
been made for departure from perfect gas. Lines of constant 
volume, cubic feet per pound mol (labeled V), and lines of 
constant pressure, pounds per square inch absolute (labeled ■P) i _ 
are plotted on S-T coordinates, with temperature scale from 

360° R (the base temperature) to 2520° R (2060° P). Zer o 

value has been assigned to entropy at the base temperature 
and 1 atmosphere (14.696 psia). 

The value at zero pressure (i.e. f for per f ect-gas condi- 
tions) of the enthalpy, called H° (British thermal units per 
pound mol) can be read from the first vertical scale to the 
right of the T-scale. This scale consists of a number of 
diagonal curved lines at the bottom of which appear values 
of P, the ratio of fuel used to that required for stoicho- 
metrically complete combustion. Por air, P = 0, the left 
edge of the H-scale is read. To allow for the effect of 
pressure on the enthalpy of air, refer to figure 2A or 2B, 
the line marked P s? 0, where the value (H - H°)/P is 
found as a function of temperature. Por example, when vol- 
ume and pressure - are 16 cubic feet per pound mol and 500 psia, 
the entropy is found to be -2.05, temperature is 740° R, 

H 0 ' = 2660 Btu per pound mol. Prom figure 2A, at 740° R and 

P = 0, (H - H b )/P = -0.085, from which H - H° * - 500 

x 0.085 — -42.5, and H = 2660 - 42 = 2618 Btu per pound mol. 
It w311 be noted that the pressure correction to H° to ob- 
tain H becomes negligible at pressures below 1 atmosphere 
or temperatures above 1000° R, and that, for many calculations, 
figures 2A and 2B may be ignored. Values of H° are based 

on H = 0 at the base temperature and 1 atmosphere.! 

Similar to H°, E° (the value of internal energy at 
zero pressure) can be read from the extreme right scale, and 
the pressure correction from figpre 3A or 3B, using the line 
P - 0. Por the previous example, when T = 740° R, 

E° = 1180 Btu per pound mol, (S - E°)/P = -0.10, 2 = 1180 

- 0 .10 x 500 = 1150. The base for 2 is consis t e nt with 

1 This gives H° a value of 7 Btu per pound mol at the 
base temperature. 
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-144PV 

that for H; 5 =? at the base temperature and 1 at- 
mosphere. 777,75 

Used as a Mixture Chart 

The use of the chart in determining the properties of an 
air— fuel mixture will next he described, without inclusion of 
the thermodynamic basis for the procedure* Suppose the mix- 
ture contains 20 percent more fuel than the stc ichiometr ic 
mixture. Then P = 1,20. The mixture may, in addition, con- 
sist in part of unpurged combustion products from the previ- 
ous cycle in an Otto cycle engine, together with some fresh 
fuel— air mixture. The weight fraction of the former in the 
total mixture is referred to as f, the ’’unpurged fraction, r 
or weight— fract ion of the fuel in the burned state. for 
Otto— cycle calculations, f may vary from 0.01 to 0.05; for 
combustion in a gas turbine cycle, f is C on entry and 
about 1,0 cn leaving. As an example, consider a mixture with 
P — 1.2 and f — 0.05, at a temperature pf 9500 R and pres- 
sure of 100 psia. Prom figure 5, one pound of air produces 
0.03545 mols of such a mixture. This is the conversion 
factor in going from the present chart to the burned charts, 
(See reference 1.) Because in many problems the oxygen of 
the air is partially or wholly combined with fuel, the ordinate 
of figure 5 is labeled ’’pound mols of mixture per pound of air 
constituents. 1 ’ Prom the main chart, at T = 950° R and 
P = 100, V a i r = 102, This volume applies for air. Since 
the p— v— T relations for air and for the mixture are somewhat 
different, a correction is read from figure 7; “ V air 

= 0.07, = 102.07, Since the uncorrected volume can be 

read tc no better than l/2 perpent , the correction at this 
temperature and pressure is obviously not worth making. 

To obtain enthalpy, first determine H° from the right- 
hand scale by moving hori?ontally along the T ^ 950— line to _ 
P - 1.2; and interpolate between the dashed— line scales cor- 
responding to an f of 0,05, finding H° = 4590. Similarly, 

B 0 = 2690 (linear interpolation and extrapolation are valid 
for obtaining values at f other than 0 or 0.C5), Prom 
figure 2A , (H - Ho)/P = -0*08 and ( K - E° ) /F - -0.1?. 

Then, H = 4590 - 0.08 x 100 4582, and t * 2690 - 0,12 X 

100 = 2678. The corrections, for this example, were barely 
worth making, 

The only additional property of interest is entropy, 
Absolute values of entropy for mixtures do not appear on the 
chart; but a change in entropy along a specified path or the 
change in other properties along a path of constant entropy 
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may be evaluated, The variation, with temperature and prei- 
sure, ef the entropy of a mixture differs from that of pure 

air because of differences in specific heats and in p-v-T 

relations for the two. Allowance for the specific-heat ef- 
fect appears in families of curves running almost vertically 
along the left side of the chart, labeled (SjqIx • $air)* 
and plotted against temperature. Allowance for the effect 
of differences in p-v-rT relations of air and of gas mixture^ 
is much less important, and may be neglected for many uses of 
the chart, (See later example.) 

As an example, suppose the last given mixture is to be 
compressed i sentropically from its temperature and pressure 
of C5Q° R and 100 paia, respectively, to a new pressure of 
600 psia, The entropy of air at the first condition is found 
to be 2.99, and the difference between the entropy of the 
mixture of interest and that of air (exclusive of a constant 
whish is independent of temperature and pressure) is read 
from the family of almost vertical lines at the left $f the 
plot. It is found to be 0.63. The sum of these is 3.62, 

Nov move upward from 950° R and 100 psia along such a path 
that the new value of the sum of S a £ r and (S m ix - S a j r ) 

is the same. This may be done by trial and-error, choosing 

various temneratures until the condition is satisfied. Mu©h 
simpler is the graphical procedure of laying the edge of a 
sheet of paper horizontally along the 950° isotherm and mark- 
ing on it the distance between the points where the line 
P =s 100 on the right, and F = 1.2, f * 0.05 on the left, 
cut the isotherm, Then, move the horizontal edge of tho 
paper vertically upward keeping the left-hand point on the 
line labeled '*F * 1.2j f * 0.05 11 until the right-hand point 

coincides with the line, P = 600, From this procedure it if 

apparent that the families of lines at the left may be thought 
of as direction lines for moving along a con st ant- entr epy 
path with a given mixture of fuel, air, and combustion prod- 
ucts. At P s* 600 it is found that T » 1435° R end 
^air =* 36. Strictly, a pressure correction to entropy should 
be made, using peripheral figure 8, The correction io appre- 
ciable only for the lower left region of the main chart (low 
temperature, high pressure). Furthermore, numerical values 
read from figure 8 indicate that, along paths of approximately 
constant entropy the pressure correction is substantially con- 
stant (within a AS of 0.005, i.e., one-tenth of a 8mg.ll 
square) and may therefore be neglected. The only remaining 
correction is to convert from V aiT . * 0 ^"mixture* the 

initial temperature the difference of the two was negligible. 
(See previous example.) At the final temperature of 1435° R 
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the correction is even less important. The work of compres- 
sion in a flow system is AH. By rounding off the result 
from the previous example, H x = 4580; R 3 - 8790, 

(H a - H° 2 )/P - 0.05; H a - H° 3 = 30, H 2 =8830. Then, 

work =8820 - 4580 = 4240 Btu. This number Is established 
to about. 10 Btu by careful use of the chart. Omission of 
the pressure corrections to H would have given the answer 
4200 Btu - that is, 1 percent lower. Treatment of the whole 
problem as though the gas were air rather than a partially 
burned fuel-air mixture would have given a final temperature 
of 1540° R and a work of compression of 4390 Btu (3.5 percent 
high) . 

CONVERSION FROM PRESENT CHART TO H BURNED MIXTURE «* CHART S 

* OF REFERENCE 1 


The present chart will be referred to as the "mpdified— 
air" chart for brevity. Given a system located on the modi- 
fied-air chart, unreacted, and of such composition that com- 
bustion would raise the temperature to a point where dissoci- 
ation is important - that iB, to a temperature not -included 
on the chart - what is the procedure? 

Cas e I Combustion occurs at constant volume. It is 
necessary to evaluate E of the burned system above the 1 older 
burned— chart base — that is, above a state in which zero value 
of internal energy is assigned to C0 2 , H 2 0(v), N 2 , 0 2 at 
100° F. From the ,, modif ied— air” chart, read E at state in 
question before reaction, and read E° at ICO 0 F (=560° R), 
To help the reader, a heavy line appears in the E— scale at 
560° R. From the same chart (fig, 5), read y = pound mols of 
mixture per pound of air constituents. Then (E — E° 56 o) y 

= sensible internal energy per pound of original air, above 
560° base, of the gas, which contains both fresh fuel— air 

and unpurged gas. The internal energies of combustion of 

these, at the base temperature, must be added to obtain E. 

■^burned chart = y ® ° seo 3 modif led— air chart - 

+ f 


+ (1 - f ) ( AE ) 


C, 560 
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The term — AE C , 5S0 is the internal energy of combustion of 

any fresh fuel in the mixture, and equals +2,151,000 Btu per 
pound mol of octene, or 

2151000 (0.01749 F) 

- = 1299 F 

28.95 

per pound air constituents?, where (0,01749 5*) is mols fuel 
per mol air, (See fig. 4*) And Q, v is the internal energy 
of combustion, at 100° F , of any unburned fuel in the unpurged 
products of a previous cycle. It Is read from the right-~hand 
side of the burned charts, or from the following table; 


00 

ft 

o 

0,9 

1.0 

1.1 

1.2 

1.5 

'z 

2 

2 

169 

336 

832 


Obviously, Q, v has no importance for lean mixtures. 
Recapitulating, • 

^burned chart = Y CE — E ° 5eo^ modif led— air chart 

+ ( 1 - f >1299 F + fQ v ( 1 ) . 

. . ✓ 

In addition to E, volume V must be known to locate 
the system on the burned chart. On the modified— air chart 
V is per mol of mixture; on the burned chart is is per pound 
of air- constituents. Figure 5 gives the conversion factor . y. 

y ^modif ied— a ir chart ^ = ^burned chart 

Cas e II. — Combustion occurs at constant pressure. The 
condition to satisfy is that H is constant. On the burned 
charts H g is the AH, to the temperature in question from 

that temperature where both E s = «r(PV) s= —NET and 

E s = ^i3MC v -( T — 560 ) f Solutions of these eauations from the 
data used in calculating the burned charts indicate that, over 
the burned— chart range of F = 0.8 to 1,5, T - 405 Is an 
adequate average base temperature, constant to within 1° F. 




IT A CA TN No* 1026 


8 


From the modified — air chart* read H at state in ques- 
tion before combustion; and read H° 405 ° ^ (from the solid 

line on the H— scale); and from figure 5, read y = pound mols 
mixture per pound of air constituents. Then (H — Ho 405 ) y 

- sensible enthalpy, per pound of original air, above 405° R 
base, cf the gases containing both fresh fuel and unpurged 
gas . 

^burned chart ^ ^ ^ 405 ^modif led— air chart 

+ ( 1 - f )( 1299 F) + f Q v ( p) 

It is to be noted that the chemical contributions to the 
value of H buraed chart involve AE C and Q v , and not 
AH C and Q,p . This is a consequence of the fact that the 
burned charts define H as B + PY , and H s as E e + PV , 
where E B * S -* . Therefore, H » H g + Q v and not 

H S + V 


Examples of Cycle Calculations 

I. Turbocomuressor -power -plant . - Consider a turbine and 
compressor, direct-coupled. The compressor takes air at 0° F 
and 8 psia, compresses it sevenfold in pressure, with an 
efficiency of 75 percent. Liquid fuel at 40° F is then 
sprayed in to give an air-fuel weight ratio of 73:1, combus- 
tion occurs at constant pressure, and the gases are expanded 
through a turbine with an efficiency of 90 percent to a pres- 
sure of 10 psia. Per pound of fuel burned per hour, what is 
the available shaft horsepower? 

At 460° R and 8 psia, find S = 2.92, E 1 ° * 696, 

(H 1 - H 1 °)/P = -0.27, H x * 696 - 0.27 X 8 - 694. 

Perfect compression to 8 x 7 = 56 psi would give 
T - 800, H = 3085 - 0.07 X 56 = 3081. Then AH would be 
3081 - 694 = 2387. Instead, because of imperfect compression 
process, 

Ah = 2387/&. 75 * 3183 

Then H 2 * 694 + 3183 =? 3877, T s = about 912° H, ( H s - H s °) 
= -0.04 x 56 = -2, H s ° = 3879, and T a = 912° H. 
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Now, add liquid fuel. From figure 4, when pounds fuel 
per pound air = l/73 = 0.01370, 


F 


0.0137 

0.06775 


0.2021 


H of liquid, per mol, from figure 9, a -13,200. and AH 
of combustion at base temperature of 360° R = 2,150,000 Btu 
per pound mol . Then 

Total H of liquid a 2,150,000 - 13,200 = 2,137,000 
From equation on figure 4, 


Hols fuel per mol air * 0.01749 X 0.2021 * 0.003535 
H s of mixture per mol of fresh mixture 


3877 + 0.003535 X 2137000 _ ^ 3g0 
1.003 54 

Mols burned mixture/mol fresh mixture is read from figure 6, 
line marked f = 0, at F = 0.2021; it is 1.011. 

When combustion occurs, H does not change. But, since 
the number of mols changed from 1 to 1.011 due to combustion, 
and since the chart basis is 1 mol 


h 3 


Ha 

1 .011 


11,265 Btu per mol burned mixture 


From chart at H 3 = 11,265, T 3 is about 1865° R, (H - H°)/P 
* +0,04 (fig. 2B), H 3 ° a 11,265 - 0.04 X 56 = 11,263, and 
T 3 a 1863° R (1403° j?) 

S 3 , air, from chart at p 3 * 56 psi and T 3 a 1863° R, 
is 9.23; Suiix - S a i r is read from the group of nearly ver- 
tical lines at the left; at F a 0.202, f = 1, and 
T 3 a 1863°, it is 0.18. The pressure correction on S is 
negligible. S 3 * 9,23 + 0,18 » 9 f 41 (plus a constant of no 
interest here ) . 
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Let the gas expand through the turbine to 10 psia. If 
perfeet expansion occurs, S(total) is still 9.41. To find 
the condition at end of this expansion, lay the edge of a 
paper hor i z on tally at T 3 * 1863 from the initial coordi nat e , 

® 3 ,alr * 9,23, to the extreme left to line 3P * 0.2, f =* 1 . 
The distance between marks on the paper now represents 
Stotal* Move the paper downward with edge kept horizontal, 

and left mark kept on the line I 1 * 0.2, f » 1 , until the 
right-hand mark intersects the line p * 10. At this point, 

T 4 '» 1205° R, H 4 * 6100. (Actually, the left-hand S-guide 

line is so nearly vertical that an approximate method is fre- 
quently better than the exact graphical one describe^). If 
the expansion path had been assumed to occur along a line of 
constant S a ^ r , a preliminary value of T 4 of 1195° R would 
be read, from which - ®air 9 0.11. a considerable er- 

ror in T 4 does not affect this quantity appreciably. Con- 
sequently, without further trial and error, S 4ja ^ r ^ 9.42 

- 0.11 « 9.30, and T 4 * 1206, H 4 ° - 6100, E 4 6100 

- 0.02 X 10 = 6100. 

Then -AH during perfeet expansion would be 11,263 

- 6100 * 5163, Actually, because of turbine inefficiency, 

-AH =* 6163 x 0.9 * 4647 
H 4 * 11,263 - 4647 * 6616 

The correction from H 4 to H 4 ° is negligible; then 

T 4 =» 1275° R 

Basis of calculation of output; 1 mol air 
Mol b fuel * 0,00354 

Pounds fuel » 0.00354 X 112.1 = 0.397 

Hols unburned mixture » 1.0035 

Mols burned mixture » 1.011 x 1.0035 

= 1.0145 

Work output of turbine » 4647 x 1.0145 * 4714 
Work into compressor * 3183 

1531 


Pounds fuel/hp-hr * 5 .. — — — = 0.66 

1631 
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The low value of this figure is a conseqixer.ce of assuming 
high values for the separate efficiencies of the compressor 
and turbine* without introducing compensating allowance for 
pressure drop in ducts and the combustion chamber. 

There may be occasion to ubo the velocity of the exhaust 
gas to obtain additional thrust. If the ex-haust gas expands 
reversibly from 10 psia to atmospheric pressure of S psia, 
the temperature changes from 1275° to 1202° R f and H° from 
6616 to 6080. Corrections to H are negligible. Then 
-AH » 636 per mol mixture. 

rs , . , , . , pounds fuel + pounds air 

Pounds mixture per mol mixture * — — — — ■ ■ — £ 

mol mixture 


0.397 + 38.95 
1.0145 


28.95 


As expected with such a lean mixture, the molecular weight 
of the products is the same as that of air. 

AH per pound mixture » » 18.52 Btu 

28.95 


This eauals 


u 

2gJ 


U * 64.4 x 778 x 18.52 * 962 feet per second 

Pounds mixture per pound fuel * 74 


Then 


Pounds thrust per (lb fuel/hr) * — 

6 

■ - - ZUl »?». = 0.66 

3600 X 52.2 


Otto-Cycle .Engine with Supercharger 

Supercharged engine operating in air at P * 10 psia 
and T = 20° F (T * 480° R), supercharging with 70 percent 
efficiency to 30 psia, after which liquid fuel is sprayed in 
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at 60° P with 20 percent excess fuel. Ungin* compression 
ratio is 6. 


Pir 8t determine H of unturned aif-fuel mixture enter- 
ing manifold. 

Isentropic compression of 1 mol of air from 10 to 30 
pounds, starting at 480° R. 

H hefore * 836 + °- 26 x 10 * 838; S 58 2 - 77 

Compression at constant - S to 30 psi gives 


T af t er 

H after 
Isentropic work * 
Actual work input 


* 667 

* 2080 + 0.12 x 30 
2084 - 838 * 1246 
a 1246/0.7 * 1780 


2084 


If no thermal losses are assumed, all the extra work goes 
into temperature in the gas, and 

H. after compression a 1780 + 838 * 2618 

from which T is approximated as 740 R 

(H - H°)/P = 0.09 

H° * 2618 - 0.09 X 30 = 2615 

T may now be determined exactly; it is 734° R 

Volume » 263 cubic feet 

At this point, inject liquid fuel at 60° P (520° R). Accord- 
ing to figure 9, 

Mol s of liquid * 1.3 X 0.01749 * 0.0210 

H of liquid a -12,200 x 0.0210 » -256 

H of mixture, per mol mixture = 2618 - 256 _ a 2313 

1.0210 

(Note that basis has changed here.) 

H° of mixture = 2313 - 0.24 x 30 * 2306, and corresponding 
temperature 690° R and V » 246, Correction on V is ^.2 

( negligible ) . 
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Now, assume that temperature in cylinder at end ex- 
haust stroke and after expansion to 10 ps i is 2100 R. 

(This is subject to check by completion of the cycle calcula- 
tions.) One mol products (f = l) at 2100° B. and 10 psi, has 

E = 9720, V ^ 2260. For the process of mixing fresh mixture 
at 30 psi through a throttling valve into combustion products 
initially at 10 psi and rising ultimately to 30 psi, while 
the volume remains constant, the enthalpy of the entering 
charge plus the initial internal energy of the combustion 

products must equal the internal energy of the final mixture. 

The basis of calculation will be 1 mol of unpurged mixture 
from previous cycle, left in engine at top dead center. 


Let x ss mols fresh charge entering at top dead center 
by throttling through intake valve and mining with unpurged 
gas, but not yet including any fresh charge entering during 
intake stroke. One mol of fresh mixture of composition 
F = 1.2 produces 1.10 mols products (fig, 6). To find the 
weight fraction of burned material in the resultant mixture, 
convert burned and unburned parts to a common state, either 
burned or unburned, and then determine the mol or volume 
fraction. By converting the fresh mixture to the burned 
state. 


1 + 1 . 1 x 

The prime indicates that this f is an intermediate value 
and not the one corresponding to the introduction of all the 
fresh charge; that comes in a later step. 


and 


■“'per mol 


per 


of 


mol 


2313x+ 9720 


resultant mixture “ 


of resultant mixture 


x + 1 

2260 
1 + x 


Given V, f 1 , and E in terms of unknown x and given 
P = 30, use the chart to solve for x by trial and error. 
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Assumed 
values of 

X 

ty _ 2260 

f i - ^ 

Values of 3 
corresp. to 7 
and f ‘ , and 

at F ' =30 

3 fr om 
above 
equa- 
tion 

Diff . 
in 
3 > 8 

1 + X 

1 + 1 , 1 X 

5 

377- 

0.154 

3320 - 3 

3550- 

-23 0 

4.6 

404- 

.165 

3800 - 3 

3638 

+ 160 

4.76 

392+ 

.160 

L 

(3590 + 10)-3 

3600 

-7 


Therefore, x * 4.76, f f » 0.160, 7 * 392, 3 * 3600, from 

which H = 5770, 

This mixture now mixes further at a pressure assumed 
constant at 30 psi, with more fresh mixture entering as the 
piston recedes to six times the previous volume. The enthalpy 
of the resultant mixture then equals the sum of the H's of 
the two components. 


6 x 2260 = 13,560 

total volume of new mixture containing \ mol of turned gases, 
4.76 mol s of fresh charge, and x x additional mole of fresh 
charge now entering. By this process, point (l) on the p-v 
diagram of the Cycle is reached. (See fig, 10.) The result- 
ant enthalpy Hj is given by 

H 1 (l + 4.76 + x x ) - 5770(1 + 4.76) + 2313 x x 

13560 

1 1 + 4.76 + xj. 


1 + (4.76 + x t )l.l 

As "before, x x is obtained "by trial and error* 


Assumed 

Y w 13560 


H corresp. 
t o V and 3 
at P=3 0 

K from 
above 
equa t i o n 

Plff. 
i n 
H ' s 

A " 5.76+Xi 


45 

267 

0.01794 

2910-6 

2704 

+ 200 

48 

252 

.0169 

2620-6 

2684 

-70 

47.3 

265 . 5 

.01716 

2690-6 

2688 

-4 


Therefore, x x » 47.3, Y x » 255,5, f * 0.0172, and 
Hi - 2688. The corresponding » 1280. 
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The next step is compression to state 2, assumed isen- 
tropic. The volume correction on V x is -0.16 (negligible). 

The entropy" 5 X is 3.37 and the entropy correction is 

-0.00013 x 30 (negligible). The isentropic compression is 
followed by measuring the horizontal distance from point 
P 1( Vj to the entr opy-dir ect i on line for F f= 1.2, f = 0.017, 

then moving the measuring device upward with left end on the 

entropy-dir ecti on line until the right end is on the line 
V 2 = 255.3/6 = 42.6. The correction on volume- is negligible. 

(As a check on the graphical procedure, 

s i + < s mix - Sair^ = 3 - 37 + °- 38 “ 3 - 73 

and 

S 2 + <S'mix - Sair) 2 = 3.74 + 1.01 = 3.75) 

From figure 8, the correction on S x is negligible. At the 
new point, P 2 * 321, T 2 - 1260, and E 2 * 4735 - 0,08 

x 321 = 4710. (Fote that the E-correction is the only one 

not negligible so far.) 

The next step is a constant-volume combustion, assumed 
adiabatic and to thermodynamic equilibrium. The temperature 
reached is so high as to make the assumption of frozen chem- 
ical equilibrium on which the present chart is based, no 
longer v.ali'd. On the other hand, the new high temperature 
makes allowance for gas imperfection unimportant. The set of 
thermodynamic-property charts presented in a previous report 
must be used at this point. As previously explained, the 
burned-mixture charts are of necessity on a different basis, 
since the number of mol s associated with a given system 
changes along various thermodynamic paths and the basis must 
be one of mass rather than number of mols, 


The cons tant— volume combustion process is one at constant 
E. From eauation (l), 

E = ( E — E° ) X y 

burned chart seo modified air chart 

+ (1 “ f X-AE c>560 ) + f 

= (4710 - 390)0.0353 + 0.983x(l299 x 1.?-) + 0,017 x 336 
- 152 +■ 1532 + 6 - 1690 Btu, E 

O 

V 3 ... = 48.6 y = 48. 6 * 0.0353 = 1,504 cubic feet 


per pound air 
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Prom the burned chart for F * 1.2, T 3 , P 3 , and S g 
corresponding to E 3 a/Lt$ V s arc 6230° R (4770® P), 

1485 pei, and 0.518. Cl% i? to be understood that the cyclf 
has been calculated assuming compression apd combustion t# 

>e adiabatic, and that the pressure and temperature are con- 
sequently high,) 

The ^ext step is isentropic expansion to V 4 * 6 x 1.504 
» 9,02. Reading from the burned chart, P 4 * 160, S 4 * 1085. 
As the exhaust valve opens, the mixture at state 4 expands 
irreversibly through the valve, but that part remaining in 
the engine expands i s ent r op i cally , doing work on the effluent 
gasee, Expansion is to P 4 = 10 pei, and from the chart 
T £ * 1900° R and V 4 f = 79. The fraction ef the mixture 
staying in the cylinder is then V 4 /V 4 * = 9,02/79 « 0.1,142. 

A 8 the piston returns to top dead center, it expels five- 
sixths of this burned residue, leaving 0.1142/6 = 0.0190 
whioh is the f, or unpurged fraction, for the cycle. This 
does not agree exactly with the value 0.0173 previously de- 
termined in steps 6-6 1 -! because the temperature of the clearr 
ance gases was assumed to be 3100° R against the value 1900 
just determined. The cycle calculation using a clearance-gas 
temperature of 1900° R instead of 2100® R ceuid be repeated, 
but the result would be found to differ negligibly from those 
just obtained. 

The work of the cycle may now be calculated. The basis 
is 1 mol of unpurged gas left in the engine at top dead cen'v 
ter. To this were adddd 


4.76 + 47,3 a? 52.06 mols fresh charge 

Each mol of this contained 1/0,03531 or 28.32 pounds fresh air. 
The work of the cycle, exclusive of the pumping loop, is 


(a* - a 4 ) - (* B - *1) 


» (1690 - 1085) 52.06 X 33.32 - (4710 - 1280)52.06 


* 892,500 - 178,500 * 714,000 Bfcu 


the pumping loop, positive in a supercharged engine, 
» - p atm> 

« (13,560 - 3260) (30 t 10) « 41,800 Btu 


The work of 
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The work put Into the supercharger (assumed mechanically 
geared to the engine) was found to he 1780 Btu per mol of 
air. Mol s of air entering 


53.06 X 28.32 
28.95 


“ scvsfe 


Work into supercharger =s 1780 x 50.95 *= 90,600 Btu 

Net work of cycle » 714,000 + 41,800 - 90,600 = 665,200 Btu 


Finally, the mols fuel entering, per mol air, * 0.01749F 
(from fig. 9 ) ; so 


Mols fuel - 50,95 X 0.01749 X 1,3 = 1.069 
Founds fuel = 1.069 X 112.13 = 119.9 

Convent ional ly , the enthalpy of the liquid fuel at a 
standard temperature is used in the denominator to express 
efficiency. The values for AH and AE for combustion are 
given in the table below the nomenclature for the main chart. 

At 560° R, AH = 2,148,000 Btu per pound mol, from vaporized 
gasoline to water as vapor. The heat of vaporization of octene 
at 560° R is obtainable from the chart, inset figure 9, by 
taking the difference between the enthalpies of the vapor and 
the liauid. It is 6100 — (—10200), or 16,300 Btu per pound 
mol. Then, AH from liquid gasoline to combus t ion. product s 
(water as vapor) is 2,148,000 — 16,300 — 2,132,000 (1904 Btu/lb) 

Then, cycle efficiency ( engine— supercharger combined) 

665200 

119.9 X 1904 
= 29.1 percent 


The obtaining of displacement volume per horsepower-hour is 
straightforward, and will not be included here. 


Massachusetts Institute of Technology, 
Cambridge, Mass., April 1945. 
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APPENDIX 

I . FUNDAMENTAL DATA 


Constants used throughout these calculation are iden- 
tical with those used in the -previous paper. The only others 
involved are constants for the Beat t i e-Br id geman equations of 
state used in this report; see later. 


Specific Heat of Octene 


Specific heats of hydrocarbon vapors may be calculated 
by a summation of a number of Einstein functions, each rep- 
resenting a bond energy in the molecule, (See reference 2.) 

■ Dobratz (reference 3) has improved the method of Bennewitz 
and Hossner by adding allowance for molecular rotation. 

Values for the specific heat of propane have been plotted, 
from the papers of Bennewitz and Rossner, Dobratz. Fugassi 
and Rudy (reference 4), and Edmister (reference 5;, for aom*. 
parison with the experimental data of Konz and Brown (refer- 
ence 6) and of Brown (reference 7). Dobratz 1 method fits the 
experimental data more closely than the others (within 1 per- 
cent at 800° K) , and furthermore predicts a high-temperature 
variation much more in accord with expectations from theory 
than any straight-line extrapolation of the experimental 
data. Consequently, Dobratz 1 correction to the calculations 
of Bennewitz and Rossner has been accepted for the calcula- 
tion of the specific heat of octene. 

The formula for C° (p=o) given by Dobratz is 


C° s 4R + aR/2 + Ev.P, + g n .. r . J ~ JHh D * Zv.q. (3) 

~ Svj. 1 A 

where 

a number of bonds permitting free rotation, C-C or 
similar 

Vj, number of bonds of type 1 

• Q<i energy associated with each bond of type 1 
n total number of atoms in molecule 
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If ms number of 0 atoms in the molecule, this becomes, 
for normal or branched olefines, 



(6 + m)R/3 + (m - 2)P X + P 2 + 2mP 3 



3 

1 


[ ( m - 2 ) + 2mQ 3 3 ( 4 ) 


For octene, for which m = 8, 


0° = 7R + 6P X + P 2 + 16P 3 + || (6^ + + 16Q 3 ) (5) 


The Sinstein function by which P x and Q,i may be deter- 
mined is 


P or Q, = R (x z -r') 

^ (e x - l) zJ 


( 6 ) 


where 
x = hu/kT 

h Planck’s constant 

V characteristic frequency 

k gas constant per molecule 

T temperature, °K 


h/k s 4.778 x 10"* 11 degree seconds 

V - c/X * 2.99796(l/X) X 10 1 0 cm/ sec/ cm 

hv _ 1 .43 24(1/ X) 

• kT T 


Wave -number s l/X 1 , 1/X gt and l/X 3 for P-energy are 990, 

1620, and 2990, and for Q-energy are 390, 845, and 1320. A 
table of Slnstein functions may be found in Phy sik. -Chem , 

Tab . ^Landol t-B8rnst eln , 1927, p. 702, fifth ed . , first en- 
larged volume . 
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o 

By applying these values to the equation, the of. ■ 

octene at 210° 23G° 250° . 990° K is determined. Then 


H ° onno v = 20.2 C 

above 200 K p 


and 


£ p - 20 2 ( C°/T ) 


Finally, E°' ( = change in internal energy above zero— base 


for E at 200° K) = 


C dT « H° - R( T - T ') 
v o 


3CC 


Entropies and Internal Energies of Other Gases . .. 

The previous report gives values of ,33 and S above 
3C0 K fcr the various "gases involved (tables II and III of 
that report). For the present purpose it. is necessary only . 
to determine values from 300° down to 200°, and pat all val- 
ues on a new base corresponding to O— value at 200° K for all 
gases. The results are given in tables I and II. 


CRITICAL CONSTANTS 


For n— octane, the critical temperature, T c , is 296° C 

'( 569° K) ; the critical pressure, P , is 24.6 atm. (See ref— 

\ ^ 

erence 8.) From ethyl ene— ethane and pr opy 1 en e— pr opan e 
cr it ical— cons tant relations (temperature down 22.4° and 5;1° 
C, respectively; pressure up 1.9 and 3.0 atm, respectively), 
assume that octene has T and which are 5° C lower and 

W C 

2 atm higher, respectively, than octane. Then ... 


T 

c , o ct ene 


564° K ; 


P 

c 


oct ene 


26.6 atm 


Considering the fact that gasoline is a mixture, the above 
procedure for determining critical constants is entirely ade- 
quate for estimating corrections due to gas imperfections, 
which in any case are small. 

For water, T c - 374.0° C(64’ 7 .2° K) 


and 


217.'? atm 
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SUMMARY OP B2ATTIB-BRJDG1MAN CONSTANTS POR &ASSS 


Ga s 

Oct ene 1 

(Air-N s 
(reference 9) 

C0 2 

(reference 10) 

H 2 0 

(reference ll) 

B P 
o c 

RT C 

■^o^c 

c 

c p c 
R T c 4 

0.1620 

.4253 
. 0667 



0.1384 

.3140 

.1438 

mm 

0.2819 

0,04070 

0.10476 

0.03377 

1 

34 . 248 

1.0763 

5.0065 

4.0686 


7 

4 

4 

4 

mam 

2.082X10 

12 x 10 

66 X 10 

951 x 10 


^Prom a private communication "by J. A. Beattie. 


For octene,. the bottom half of the table was obtained 
from the top half and from ths previously given critical data, 
as follows: 


0 



0,08206 X 564 

= 0.2819 liters/gm mol 

26.6 / 

(0 . 08206 x 5 64 ) S 

■ = 34 . 248 (liter s/gm mol) 2 (atm) 

26.6 

0.08206 X 564 4 _ , 

= 2.082 X 10 (liters/gm molK *0 

26.6 


The Beat t i e-Br i dgeman constants for water vapor were obtained 
in like manner. 

It was assumed that all diatomic species (CO, 0 S , H s , OH, 

NO), because of the minute quantity present, had the same B.B. 
constants as Air-^N g . 
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II. CALCULATION OF TEKPSRATURS AT WHICH H = 0 
ON BURNED-MIXTURE CHARTS 

F = 0.8 Composition of gas is from page 55 of M,I.T. file- 
notes. 


Composition (AE,200 to 300° K)/l00 *= MC V N.HC V 

(cal/ gm mol ) 


Ns 

1 

4.95 

4.9500 

NO 

.00119 

5.21 

. 0062 

o s 

.05248 

5.02 

,2635 

COg 

.14176 

6 .40 

.9080 

HgO 

.14154 

6.00 

.8500 

CO 

.00001 

4,97 

_ — — 

H s 

— 

4.78 

--- 

OH 

.00045 

5.21 

.0027 


EN=1. 33743 


2NMC V =6 .9804 


- 

-NET = 2NMC V (T _ T 0 ) 


-1 

.33743 x 1.987 

X T = 6 .9804 (T - 560) = 

6.9804T - 3910 


T 

3910 Q 



* = 406° 

E 



6.9804 + 2.655 


F = 

1^5 




Composition 

KC V 

N.MC V 

N e 

1.0 

4,95 

4.9500 

C 0 2 

.09667 

6.40 

. 6185 

h 3 0 

.16897 

6.00 

1,0140 

CO 

.16900 

4.97 

.8400 

S 3 

. 09669 

4.78 

.4625 

H 

. 00002 

3,00 

.0001 


2N=1. 53134 


ZNMC V =7 .8851 


-1.53134 

X 1 .987 X T = 7 , 8851 (T 

- 560) 



-3 .04 T = 7.8851 T ~ 

4415 



_ 4415 o 




T « * 404 P 




10,925 
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Spread of but 2° V between t£ie two extremes of composition 
suggests an average T = 405° R at which H = 0. ThiB pro- 
duces a maximum error in H of about 7 Btu per mol, or about 

0.25 Btu per pound, whieh on the burned charts is one- 
thirtieth of one small division on the H-scale. A heavy line 
has befen drawn on the Upscale of the Modified Air chart at 
T = 405 R. 
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E 


w 




o 

W 

5 a 

00 

SO 

OH 

0 ,H 

0 

0 

0 

0 

0 

0 

299 

235 

248 

263 

262 

149 

600 

478 

497 

521 

521 

298 

1221 

97^ 

997 

1047 

1032 

596 

2508 

1975 

2025 

2101 

2046 

1192 

3907 

2986 

311 8 

3237 

3067 

1788 

5*428 

4017 

1(279 

IA39 

4 n 6 

2384 

7071 

5083 

5495 

5691 

5198 

2980 

8816 

6188 

6751 

6979 

6338 

3576 
















TABLE II.- CONSTANT -PEES SUBE CHANG* OP iKXEOFX WITH IW-3WJUBK , S° 
(oal/gm mol °K or Btu/lb mol °R) 

[Zero value assigned at 200 °K] 


I 

(°K) 

Air 

Octane 

Air-Jig 

o 2 

u ■■ . 

COg 

H a 0 

...... 

Hg 

CO 

NO 

OH 

0,H 

r 

200 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

250 

1.544 

6.176 

1,542 

1.553 

1.809 

1.773 

1.495 

1.553 

I.615 

1.609 

1.11 

1 

300 

2.815 

12.141 

2.812 

2.826 

3-3*4 

3.230 

2,742 

2.821 

2.921 

2.913 

2.02 

4oo 

4.824 

23.905 

4.812 

4.869 

6.090 

5.556 

4.739 

4.S29 

4-978 

4.98 

3-45 

0 

0 

7- 714 

46.987 

7.671 

7.876 

10.3S9 

0.959 

7.573 

7.717 

7-937 

7.84 

5.46 

! soo 

9.S54 

68.579 

9-778 

I0.l4i 

13.790 

11.531 

9.594 

9.859 

IO.I36 

9.88 

6.89 

1000 

11-585 

88.277 

11.482 

11.972 

16.617 

13.657 

11.190 

11.596 

11.921 

11.49 

8.00 

1200 

13.049 


12.927 

13.510 

19.034 

15,508 

12.523 

13.065 

13-39 

12.84 

8.90 

! 1400 

•*. — 

14.319 


14.181 

14.837 

21.143 

17.161 

13.680 

14.339 

14.69 

14.01 

9.67 


CO 

o» 
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